Quantum correlations dynamics under different non-markovian environmental models 
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We investigate the roles of different environmental models on quantum correlation dynamics of 
two-qubit composite system interacting with two independent environments. The most common en- 
vironmental models (the single- Lorentzian model, the squared-Lorentzian model, the two-Lorentzian 
model and band-gap model) are analyzed. First, we note that for the weak coupling regime, the 
monotonous decay speed of the quantum correlation is mainly determined by the spectral density 
functions of these different environments. Then, by considering the strong coupling regime we find 
that, contrary to what is stated in the weak coupling regime, the dynamics of quantum correlation 
depends on the non-Markovianity of the environmental models, and is independent of the environ- 
mental spectrum density functions. 

PACS numbers: 03.67.Mn, 03.65.Yz, 42.50.-p, 71.55.Jv 



I. INTRODUCTION 



Until recently a lot of interest has been devoted to 
the definition and understanding the quantum aspects of 
correlation in a composite system. The discovery that 
mixed separable (unentangled) states can have nonclas- 
sical correlation [1-4] and such states provide computa- 
tional speedup compared to classical states in some quan- 
tum computation models [4,5]. In order to quantify the 
quantumness of the correlation contained in a bipartite 
quantum state OUiver and Zurek [3] proposed a measure 
for quantum correlation known as quantum discord (QD) 
and based on a distinction between quantum information 
theory and classical information theory. A recent result 
that almost all quantum states have a nonvanishing QD 
[6] shows up the relavance of studying such correlation. 

Besides the quantification of quantum correlations, an- 
other important problem is the behavior of these correla- 
tions under the action of decohcrence. The phenomenon, 
caused by the injection of noise into the system and aris- 
ing from its inevitable interaction with the surrounding 
environment, is responsible for the loss of quantum co- 
herence initially present in the system. Recently it was 
noted that the dynamical behaviors of QD in the pres- 
ence of the Markovian [7,8] decoherence decay exponen- 
tially in time and vanish only asymptotically [9,10], con- 
trary to the entanglement dynamics where sudden death 
may occur [11-18]. In these above studies, the quantum- 
ness of correlation is more robust to the action of the 
environment than the entanglement itself. In particu- 
lar, Refs. [19,20] have discovered that the QD can be 
completely unaffected by Markovian depolarizing chan- 
nels or non-Markovian depolarizing channels for long in- 
tervals of time, and this phenomenon has been observed 



experimentally [21]. As Refs. [19-21], it is of interest 
to find a certain environmental model that the quantum 
correlation can be unaffected by decoherence as much as 
possible. 

In this article, we will concentrate on the question: 
what kind of local environmental model can make the 
initial quantum correlation more robust in the dynamics 
process? Wc consider a noninteracting two-qubit system 
under the influence of two independent environments. 
The most common environmental models (the single- 
Lorcntzian model, the squared-Lorentzian model, the 
two-Lorentzian model and band-gap model) are studied. 
By analytical and numerical analysis we find that, for the 
weak coupling regime, the monotonous decay speed of the 
two-qubit QD is mainly determined by the spectral den- 
sity functions of these different environments. The two- 
qubit QD in the single-Lorentzian (band-gap) environ- 
ment is more robust than in the squared-Lorentzian (two- 
Lorcntzian) environment under the resonant and near 
resonant conditions. While for the far off-resonant con- 
dition the two-qubit QD in the single-Lorentzian (band- 
gap) environment decreases much more faster than in the 
squared-Lorentzian (two-Lorentzian) environment. How- 
ever, by considering the strong coupling regime we find 
that, the two-qubit QD is more robust in the squared- 
Lorentzian (two-Lorentzian) environment than in the 
single-Lorentzian (band-gap) environment, either under 
the near resonant or far off-resonant condition. In this 
case, the dynamics of QD mainly depends on the non- 
Markovianity of the environmental models, and is in- 
dependent of the environmental spectrum density func- 
tions. 



II. THEORETICAL MODEL AND DYNAMICS 
OF TWO-QUBIT SYSTEM 
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Considering a model consisting of two qubits A and B, 
each interacting with a zero-temperature bosonic envi- 
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ronmcnt, denoted a and b, respeetively, we assume that 
each qubit-environment system is isolated and the en- 
vironments are initially in the vacuum state while two 
qubits are initially in an quantum correlated state. A 
specific system which consists of two independent two- 
level atoms interacting with an multi-mode environment 
respectively has been chosen in this paper. Since each 
atom evolves independently, we can learn how to charac- 
terize the evolution of the overall system from the atom- 
environment dynamics. The interaction between an atom 
and an N-mode environment in the rotating-wave ap- 
proximation is given by Hq + Hf^^, which, in the basis 
ilgg), Ige), leg), |ee)}, reads 

TV 

H^Q = woCT+fTi + ^ Wfcaj^afc, (1) 
fc=i 

TV 

HLt = 5If'^('^>'^ + ^-4)> (2) 

k=l 

here aj,, ak are the creation and annihilation operators 
of quanta of the environment (a or b), = |ej)(gj|, 
(T^ = \gj){ej\ and ujj are the inversion operators and 
transition frequency of the j-th atom (j=^, B and here 
ujA = wb = Wo ) ; and gt are the frequency of the mode 
k of the environment and its coupling strength with the 
atom. To illustrate the roles of the different environ- 
mental models on quantum correlation dynamics of two 
atoms, we assume that two atoms interact off-resonantly 
with their structured environment, whose spectral den- 
sity function D{uj) provides a complete characterization 
of the evolution for single-Lorentzian, two-Lorentzian, 
band-gap and squared-Lorentzian environments. 

In order to find the atom-environment dynamics, we 
solve the master equation by using the pseudomodc ap- 
proach [22,23]. This exact master equation describes the 
coherent interaction between the atom and the pseu- 
domodes in the presence of the decay of the pseudo- 
modes due to the interaction with a Markovian reser- 
voir [24]. The number of the pseudomodes relies on 
the shape of the environmemt spectral density func- 
tion, (a) For the single-Lorentzian environmental model 
D{uj) = )^+(r/2)^ ' there has only one pole in the 

lower half complex plane, the atom interacts with one 
pseudomodc which leaks into a Markovian environment. 
So the exact dynamics of the atom interacting with a 
single-Lorentzian structured environment is contained in 
the following pseudomode master equation 

dp r 

— = —i[H\ p] — — [a)ap — 2apa) + pa)a\^ (3) 

where 

= cjocr+crl + oj^a^a + n{(T^_^_a + ciia^), (4) 

with p is the density operator for the j-th atom and the 
pseudomodc of the structured reservoir, a and a' are the 



annihilation and creation operators of the pseudomode. 
The constants Wc and F are, respectively, the oscillation 
frequency and the decay rate of the pseudomode and they 
depend on the position of the pole z = Wc — «F/2. The 
j-th atom interacts coherently with the pseudomode (the 
strength of the coupling il). 

(b) According to the two-Lorentzian environmen- 
tal model, the environment spectral density func- 
tion is simply a sum of two Lorentzian functions 

^(^) = ^ H'-'-J+ir,/2Y^ + ^2 (^-^,)l;(r,/2)^ ^ ^here 
the weights of the two Lorentzians are such that Wi + 
W2 = 1. There are two poles in the lower half complex 
plane, the atom interacts with two pseudomodes (ai and 
02) which leak into a Markovian environment (Fi and F2 
are the decay rates), respectively. This time the poles 
are located at Zi = uic — iTi/2 and 22 = — «F2/2, so 
the exact master equation for the atom-environment dy- 
namics in the two-Lorentzian environmental model can 
be written 

^ = - i[H^,p]-^{a\aip-2aipa\+ palai) 

- -2-(4a2P - 2a2p4 -t- p4«2), (5) 

and here 

— ujQ (T-'^a-L + ujca\ai + ujca\a2 + il-\/VFi(cr:^ai 
-I- (Jla\) + n^/w'2{a^+a2+(Jlal). (6) 

(c) Next we give an idealized model of a band gap (or 
photon density of states gap) D{lu) ~ , ~ 

^^^^ I ' — in which both Lorentzians are centered at 

the same frequency, the second is given a negative weight, 
and the weights of the two Lorentzians are such that 
Wi — W2 ~ 1 and F2 < Fi ensure positivity of D{uj). 
There also have two poles in the lower half complex plane 
as the two-Lorentzian model, the two poles are located 
at ujc — iFi/2 and uJc — iT2/2, so there are also two pseu- 
domodes ai and 02 with deacy rates T[ = W1F2 — W2F1 
and F2 = WiTi — W2T2 respectively. The j-th atom 
does not couple to the first pseudomode oi at all, it 
only interacts coherently with the second pseudomode 
02 (the strength of the coupling fl) which is in turn 
coupled to the first one (the strength of the coupling 
V = ^/W[W^{ri — F2)/2), and both pseudomodes are 
leaking into independent Markovian environments. The 
exact pseudomode master equation associated with the 
band-gap model is given by 

dp F' 

— = - i[H^ , p] - -^[a{aip - 2aipal + pa\ai] 
F' 

- y[°2^2^' ~ 2a2p4 + P4a2], (7) 

where 

= ojQ(j-'^ai_ + Wcftjai + uJcala2 + fl{a\(ji_ 

+ a2C7^+) + V{a{a2 + aial). (8) 
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(d) The environment spectral density function of the 

squared-Lorentzian model is D{uj) = [(^_^ )2^(r/2)2]g ' 
for which we will find that there exist two pseudomodes 
fli and a2, and the j-th atom only couples to the second 
pseudomode a2 (the coupling constant il) which inter- 
acts with the first pseudomode ai (the strength of the 
coupling V = r/2). Different from the band-gap model, 
only the first pseudomode will show any decay to the 
Markovian environment with decay rate T, the second 
pseudomode which is directly coupled to the j-th atom 
does not decay in this model. So the dynamics of the j- 
th atom and two pseudomodes obey the following master 
equation 

^ = -i[H^,p\ - T[a\aip - 2aipal + palai], (9) 
with the Hamiltonian 

= ojQd-'^cr'L + ujca\ai + L0ca\a2 + Vt{a\a'^_ 

+ a2(T^_^) + V{a\a2 + aia\). (10) 

In order to analyze the roles of the different envi- 
ronmental models on quantum correlation dynamics of 
two atoms, we consider the above four environmental 
models, respectively, i.e., the single-Lorentzian model, 
two-Lorentzian model, band-gap model and squared- 
Lorentzian model. According to the above analysis, 
the spectral density functions of single-Lorentzian model 
and squared-Lorentzian model have a same parameter 
r, and the two-Lorentzian model and band-gap model 
both contain two Lorentzians, the same parameters [Wi, 
W2, Fi and appear in the spectral density func- 
tions of them. So in this paper we will mainly com- 
pare the difference in quantum correlation dynamics 
of two atoms between the single-Lorentzian model and 
squared-Lorentzian model, as well as bewteen the two- 
Lorentzian model and band-gap model. For an ini- 
tial state of the total system p{0)ABab — with 
= {cos0\gg)AB + sm9\ee)AB) ® |0)a|0)6, and here 
9 e [0, tt], |e) and \g) arc the excited state and ground 

state of atoms, \0)a.b = Ylk=i \^k)a,b is the vacuum state 
of the environment a, b. Then the evolutional density 
matrix p(t) of the total system in different environmen- 
tal models can be acquired respectively by solving the 
above master equations (from Eqs. ([3]) to ([T0|) ). Tracing 
out the pseudomode degree of freedom, wc obtain the 
reduced density matrix pAsit) of the atomic system in 
these four different environmental models. 

The measure of total quantum correlations used here 
is the quantum discord (QD) [3]. In all cases investi- 
gated in this paper the reduced density matrix for the 
atomic system pAsit) in the basis {\gg), \ge), \eg), |ee)} 
has an X structure defined by its elements pi2 — P13 = 
P24 = P34 = 0, P22 = P33 and pi4 = pIi- For this 
X class of density matrix, QD can be calculated ana- 
lytically [25]: QD{pab) = ^(pb) + E'=o log2 + 
min{5. J [S'(p|{_Bi})], where Xj is the j-th eigenvalue of 




FIG. 1: (Color online) Time evohition of the atomic QD as 
a function of the dimensionless quantity fit under the atom- 
pseudomode resonant condition (A = ojc — 1^0 = 0), with 
9 = 7r/3. (a) and (b) the single-Lorentzian and squared- 
Lorentzian as the environmental models; (c) and (d) the two- 
Lorentzian and band-gap as the environmental models. 

the density matrix pAsit). Here S{pb) denotes the von 
Neumann entropy oi ps = TtaPab and S{p\{Bi}) is the 
quantum conditional entropy with respect to a von Neu- 
mann measurement {Bi\ for subsystem B. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In Figs. 1(a) and 1(6), by considering the atom- 
pseudomode resonant condition (A = Wc — wq = 0) 
and choosing the single-Lorentzian model and squared- 
Lorentzian model as the environmental spectral density 
functions, we plot the time evolution of QD for two 
qubits as function of the dimensionless quantity 51t in the 
weak coupling regime and the strong coupling regime, 
with 6 = 7r/3. In the weak coupling regime QD decays 
only asymptotically to zero while in the strong coupling 
regime the atomic QD presents damped oscillations. A 
comparison between the dark solid curve and the red 
dashed curve in Figs. 1(a) and 1(6) reveals that for 
the weak coupling regime, the atomic QD due to the 
single-Lorentzian environmental model is more robust 
than the case of the squared-Lorentzian model, but in 
the strong coupling regime, the atomic QD attenuates 
more slowly in the case of the squared-Lorentzian model 
than the single-Lorentzian model. Through comparing 
the atomic QD dynamics in the two-Lorentzian envi- 
ronmental model and band-gap environmental model 
(as shown in Figs. 1(c) and 1(d)), we acquire that in 
the weak coupling regime, the QD can decrease to a 
long-time asymptotic value in the band-gap model while 
for the case of the two-Lorentzian model the QD can 
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reduce eventually to zero. However, according to the 
strong coupling regime, the atomic QD is more robust in 
the case of the two-Lorentzian model than the band-gap 
model. 

Then, in order to investigating the effects of different 
environmental models on the atomic QD under the 
atom-pseudomodc near resonant and far off-rcsonant 
conditions, wc analyze the evolution behavior of the 
QD in the weak coupling regime, by the comparison of 
two cases: near resonance condition (A = 0.20.) and 
far off-resonance condition (A = SfJ), with 9 = it/ 3. 
For the case of near resonance, as shown in Figs. 
2(a) and 2(c), one could find that the atomic QD in 
the single-Lorentzian (band-gap) environment is more 
robust than in the squared-Lorentzian (two-Lorcntzian) 
environment. However, an opposite result that the 
atomic QD in the single-Lorentzian (band-gap) en- 
vironment decreases much more faster than in the 
squared-Lorentzian (two-Lorcntzian) environment are 
obtained for the far off-resonant condition, clearly seen 
in Figs. 2(5) and 2(d). What is the physics behind the 
phenomena? In this part we try to give an enlightening 
discussion for this problem based on these environmental 
models. Let us review the spectrum density functions 
of these models, as shown in Fig. 3. The center part of 
the spectrum density function of the single-Lorentzian 
(band-gap) environment is much smaller than of the 
squared-Lorentzian (two-Lorentzian) environment. In 
contrast, the parts which are far from the center are 
larger in the single-Lorentzian (band-gap) environment 
than in the squared-Lorentzian (two-Lorentzian) envi- 
ronment. Thus, to determine in which environmental 
model the atomic QD is more robust in the weak 
coupling regime, we can compare the spectral density 
functions of these different environments: the decay 
behavior of the atomic QD is determined by the modes 
of the spectrum which are resonant with the atoms: the 
monotonous decay speed of the QD decreases as the 
density of these modes decreases. 

In this paper, we also understand the influences 
of different environments on the atomic QD in the 
strong coupling regimes which satisfy F = O.llfJ in the 
single-Lorentzian environment and squared-Lorentzian 
environment, and Fi = O.llil, F2 — 0.0117 in the 
two-Lorcntzian environment and band-gap environment. 
In Fig. 4, we acquire that the periodically oscillating 
decay speed of the atomic QD in the squared-Lorentzian 
(two-Lorentzian) environment is slower than in the 
single-Lorentzian (band-gap) environment, either under 
the atom-pseudomode near resonant or far off-resonant 
condition. That is to say, in the strong coupling regime 
the QD is more robust in the squared-Lorentzian (two- 
Lorentzian) environment than in the single-Lorentzian 
(band-gap) environment. In what follows, we will give 
a simple interpretation for why this finding in the 
strong regime is different from the results in the weak 
coupling regime. First, taking the spectrum density 
function D{uj) of the above four environmental models 
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FIG. 2: (Color online) Time evolution of the atomic QD as a 
function of the dimensionless quantity Qt under the atom- 
pseudomode near resonance regime and far off-resonance 
regime, with 6 = tt/3. (a) and (b) the single-Lorentzian 
and squared-Lorentzian as the environmental models, with 
F = HQ; (c) and (d) the two-Lorentzian and band-gap as the 
environmental models, with Fi = llSl, r2 = SI. 
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FIG. 3: The density of the spectrum D{iu) as a function 
of the dimensionless quantity {uj — ojc) in the weak-coupling 
regime, (a) for the single-Lorentzian and squared-Lorentzian 
as the environmental models, with F = (b) for the two- 
Lorentzian and band-gap as the environmental models, with 

Fi = im, Vi = n. 



into account in the strong regime, we note that the 
discrepancy among them is very minor, as shown in 
Fig. 5. So from the spectrum density function to give a 
construction is not feasible. 

However, according to the previous works [26-28], 
wc know that there exists the non-Markovianity of 
environment in the strong coupling regime, and the non- 
Markovian effect of environment can play an important 
role on the dynamics of the qubits system. Therefore, we 
will show the degree of the non-Markovian behavior of 
the dynamics processes in these different environmental 
models. In Rcf. [26], Breuer el al. suggest definition a 
measure N{^) for the non-Markovianity of the quan- 
tum process $(t) by means of the relation N{^) = 
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FIG. 4: (Color online) Time evolution of the atomic QD as a 
function of the dimensionless quantity Sit under the atom- 
pseudomode near resonance regime and far off-resonance 
regime, with 9 = 7r/3. (a) and (b) the single-Lorentzian 
and squared-Lorentzian as the environmental models, with 
r = O.llil; (c) and (d) the two-Lorentzian and band-gap as 
the environmental models, with Fi = O.llSl, F2 = O.OISI. 
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FIG. 5: The density of the spectrum D{uj) as a function of 
the dimensionless quantity (uj — Uc) in the strong-coupling 
regime, for the single-Lorentzian and squared-Lorentzian as 
the environmental models, with F = 0.11f2; and for the two- 
Lorentzian and band-gap as the environmental models, with 
Fi = O.im, F2 = O.Olfi. 



inaxpj,2(0) Ei[-D(pi(&j),P2(foO) - D{pi{a.i),p2{ai))]. To 
calculate this quantity one first determines the total 
growth of the trace distance over each time interval 
(ai,6i) and sums up the contributions of all intervals. 
Then N{<^) can be obtained by determining the max- 
imum over all pairs of initial states. Taking the far 
off-resonance as an example, the analytical expression 
of the trace distance in the atom-environment dynamics 
process is D^{pi,p2) = here with b{t) represents 

the amplitude damping of the excited state |e), and the 
pair of initial states pi{0) = |e)(e| and p2(0) = \g){g\ 
which optimize the total increase of D'^ipuP2). Thus, 



far off-resonanl (A^sni 

Single-Lorentzian Model 

- - Squared-Lorentzian Model 




far off-resonant {i=En) 
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FIG. 6: (Color online) Time evolution of the trace distance 
{pi, P2) as a function of the dimensionless quantity Q,t 
under the far off-resonance regime, with pi(0) = |e)(e| and 
P2(0) = (a) for the single-Lorentzian and squared- 

Lorentzian as the environmental models, with F = O.llSl; (b) 
for the two-Lorentzian and band-gap as the environmental 
models, with Fi = O.llfi, F2 = O.OIH. 



we can qualitatively and intuitively compare the non- 
Markovianity due to the different environment models 
by the time evolution of the trace distance. Fig. 6 shows 
the trace distance D^{pi,p2) as a function of fit for 
A = 8fl, namely, the far off-resonance regime. It is very 
interesting to note that the amplitudes of D^{pi,p2) 
caused by squared-Lorentzian (two-Lorentzian) envi- 
ronment are much wider than those caused by the 
single-Lorentzian (band-gap) environment. In other 
words, the non-Markovianity of the squared-Lorentzian 
(two-Lorentzian) environment is much stronger than the 
single-Lorentzian (band-gap) environment. This finding 
leads to a clear interpretation for the result obtained by 
Fig. 4: the atomic QD in the strong coupling regime 
is determined by the different degree of environmental 
non-Markovianity, and is independent of the spectrum 
density function D{u]). 

In conclusion, we have studied the quantum correla- 
tion dynamics in the different decoherence environments, 
and considered a two-atom system interacting with two 
local, independent environments, modeling several 
common noise sources: the single-Lorentzian model, the 
squared-Lorentzian model, the two-Lorentzian model 
and band-gap model. For the weak coupling regime, 
it is clear to realize that the atomic QD in the single- 
Lorentzian (band-gap) environment is more robust than 
in the squared-Lorentzian (two-Lorentzian) environment 
under the resonant and near resonant conditions. But for 
the far off-resonant condition the opposite result shows 
that the atomic QD in the single-Lorentzian (band-gap) 
environment decreases much more faster than in the 
squared-Lorentzian (two-Lorentzian) environment. How- 
ever, for the strong coupling regime, the atomic QD is 
more robust in the squared-Lorentzian (two-Lorentzian) 
environment than in the single-Lorentzian (band-gap) 
environment, cither under the atom-pseudomodc near 
resonant or far off-resonant condition. Finally we note 
that we study here only the two-atom system interacting 
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with their independent environments. An important 
future investigation will be the study of the effects of 
these different environmental models on the dynamics 
of the two-atom system under a common environment, 
where quantum correlations can be created in the 
system through nonlocal interactions mediated by the 
environment. 
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